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Abstract 

We consider a class of Hamiltonian PDEs that can be split into a linear un- 
bounded operator and a regular non linear part, and we analyze their numerical 
discretizations by symplectic methods when the initial value is small in Sobolev 
norms. The goal of this work is twofold: First we show how standard approxima- 
tion methods cannot in general avoid resonances issues, and we give numerical 
examples of pathological behavior for the midpoint rule and implicit-explicit in- 
tegrators. Such phenomena can be avoided by suitable truncations of the linear 
unbounded operator combined with classical splitting methods. We then give a 
sharp bound for the cut-off depending on the time step. Using a new normal 
form result, we show the long time preservation of the actions for such schemes 
for all values of the time step, provided the initial continuous system does not 
exhibit resonant frequencies. 
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1 Introduction 

In this work, we consider Hamiltonian partial differential equations whose Hamil- 
tonian functional H = Hq + P can be split into a linear operator Hq and a non- 
linearity P{u). Typical examples are given by Schrodinger equations and wave 
equations (see [1, 9, 2, 6, 7] for similar framework). 

The understanding of the geometric numerical approximation of such equa- 
tions over long time has recently known many progresses, see [6, 7, 3, 8]. In 
essence these results show that under the hypothesis that the initial value is 
small and the physical system does not exhibit resonant frequencies, then the 
numerical solution associated with splitting methods induced by the decompo- 
sition of H will remain small for very long time in high Sobolev norm, under a 
non resonance condition satisfied by the step size h < h^. An analysis of this 
condition ensures the absence of resonances only if a strong CFL assumption 
is made between h and the spectral parameter of space discretization. Outside 
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the CFL regime, splitting methods may always exhibit resonances, in the sense 
where for some specific values of the stepsize, the conservation properties of the 
initial system are lost. 

To avoid these problems - which would require in practice the exact knowledge 
of the whole spectrum of the linear operator - it is traditionally admitted that 
the uses of implicit and symplectic schemes should help a lot. Indeed for linear 
PDEs for which splitting methods induce resonances (see [5]), it has been recently 
shown (see [4]) that the use of implicit-explicit integrators based on a midpoint 
treatment for the unbounded part allows to avoid numerical instabilities. 

In this work, we show that in the nonlinear setting, it is in general hopeless 
to find a numerical method avoiding resonances unless the linear part of the 
equation is suitably truncated in the high frequencies. By numerical examples, 
we actually show how the midpoint rule applied to a Hamiltonian PDE either can 
yield numerical resonances and destroy the preservation properties of the initial 
system, or on the opposite break physical resonances and produce unexpected 
long time actions preservation. 

We then show that if the linear operator is truncated in high frequencies, 
splitting methods yield numerical schemes preserving the physical conservation 
properties of the initial system without bringing extra numerical resonances be- 
tween the stepsize h and the frequencies of the system. In the case where this 
high frequencies cut-off corresponds to a CFL condition, we somewhat give a 
sharp bound for the CFL condition to avoid resonances extending the results in 



Roughly speaking, the phenomenon can be explained as follows: for splitting 
schemes, resonances reflect the needle of a control of a small divisor of the form 



is a (signed) sum of the frequencies ujj of the linear operator Hq indexed by a 
countable set of indices j (the sign it depends of the index j^, see below for 
further details). The control of these small divisors up to an order r ensures the 
control of Sobolev norms of the numerical solution on times of order e"'' for an 
initial value of order e. 

Such a control can be made when hQ{j) < vr where (1.1) is essentially equiv- 
alent to hi}{j) which can be controlled using classical estimates inherited from 
the generic assumption of absence of resonances in the physical system (see for 
instance [9, 1, 2]). Such an assumption will be satisfied under a CFL condition 
bounding the frequencies in terms of the space discretization parameter, but it 



[6, 7] and [3, 8]. 



exp(i/iJ7(j)) — 1 




where j = {ji, ■ ■ ■ ,jr) is a multi-index, and where 



0(jf) = ±iOj-^ ± ■ ■ ■ ±UJj, 
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is worth noticing that it depends a priori on r which is an arbitrary parameter 
(see [6, 7, 8]). 

Outside this CFL regime, resonances can always occur when hVl{j) ~ 21tt for 
some £ E Z ^ {0}. In [7], it is shown however that such situation is exceptional, 
as it corresponds to very specific values of the time-step h. 

In the case of a general symplectic Runge-Kutta method or for implicit- 
explicit split-step method based on an implicit symplectic integrator for the linear 
part, the small divisors to be controlled take the form 

exp(i*(/i,j))-l (1.2) 

where 

and where is associated with the RK method. For the midpoint rule, the 
corresponding function is 

V'(x) = 2arctan (— ). (1-3) 

Hence for these methods the control of the small divisors depends on the rela- 
tion ^(h,j) = 2TTi which is a nonlinear relation between h and the frequencies 
u). Of course a very restrictive CFL hypothesis stating that the hcoj^, are small 
(e.g. of order h) ensures the control of these small divisors, but in general res- 
onances occur. On numerical examples, it is indeed possible to exhibit h such 
that ^{h,j) = while ^{j) 7^ (numerical resonances). Moreover, the some- 
what opposite situation is also possible: the presence of physical resonances (i.e. 
j such that r2(ji) = 0) are broken by the action of h, yielding a ^{h,j) 7^ and 
unexpected regularity preservation. We refer to the section devoted to numerical 
examples for concrete illustrations. 

To avoid resonances in general situations, it seems very difficult to avoid the 
two following restrictions: 

• 'ijj{x) is linear, i.e. the numerical schemes is base on a splitting scheme or ex- 
ponential integrator. If this is not the case, we believe that the pathological 
behaviors described above can always be observed. 

• There is a frequency cut-off in order to avoid situations where hVt[j) = 2^7r 
with ^ 7^ 0. Note that this does not exactly correspond to a CFL condition, 
as high frequencies are allowed to exist without restriction, but the action 
on the linear operator on these high frequencies is cancelled. 

In this work we give a sharp explicit bound for the cut-off in order to avoid nu- 
merical resonances and we prove a new normal form result for the corresponding 
splitting methods, yielding the preservation of the regularity of the numerical 
solution for a number of iterations of the form 

n < e-^h~^ + e-'+^ (L4) 
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where r and are given parameters, and where e denotes the size of the initial 
value in Sobolev norm. In comparison with [6, 7, 3, 8], the difference is that 
the CFL condition is only imposed on the linear operator, and moreover this 
condition is independent of the approximation parameters r and A'^, and is sharp 
in the sense where we can in general exhibit numerical resonances if it is violated. 

As a conclusion, we would like to stress the following: the action of the 
midpoint rule (or of any symplectic RK methods) on the operator Hq is equivalent 
to a smoothing in high-frequencies which amounts in some sense to a cut-off. 
What we show here is that the nonlinearity of the smoothing function (essentially 
based on the function arctan(rE)) introduces in general numerical instabilities. To 
avoid them, the user should better make the cut-off by himself! Note that this 
only requires the a priori knowledge of bounds for the growth of the eigenvalues 
of the linear operator Hq. 

2 Numerical examples 
2.1 Actions-breaker midpoint 

We consider the cubic nonlinear Schrodinger equation 

idtu = — Au + V -ku^ \u\^u 

on the torus = M/27rZ. We consider V the function with Fourier coefficients 
V{k) = — iQ^2k'^ • "^^^ frequencies of the operator Hq in Fourier basis are hence 

We first consider the implicit-explicit integrator defined as 

^,"+1 = R{-ihHo) o <^^(n") (2.2) 

where are approximations of the exact solution u at discrete times tn = nh, 
and where 

is the stability function of the midpoint rule, and where 

Lpp{u){x) = exp{—ih\u{x)\'^)u{x) 

is the exact solution of the nonlinear part. Note that the unperturbed integrator 
(i.e. without the nonlinear term) can be written in terms of Fourier coefficients 
(see [4]) 

R(—ihH())k = '- — = exp(2z arctan(/ia;fc/2)). (2-3) 

1 -|- ihujk/2 
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We consider as initial value the function 



0.1 



u 



+ 0.05(2e 



2ix 



(2.4) 



2 — 2 cos(x) 

We use a collocation space discretization with K = 100 Fourier coefficients. 
In figure 1 and 2, we plot the evolution of the sum of the numerical actions 
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r,n 1 2 



for all /c G N in logarithmic scale with respect to the time. 
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Figure 1 : Plot of the actions for non-resonant step size (mid-split integrator) . 
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Figure 2: Plot of the actions for resonant step size (mid-split integrator) 



In Figure 1, we use the stepsize h = 0.13, and we observe the long-time 
behavior of the actions. In Figure 2, we use a stepsize h = 0.1278... such that 

arctan(/ia;2/2) -|- arctan(/ia;5/2) — arctan(/itiJ_7/2) = (2.5) 

and we observe energy exchanges between the actions. This corresponds to the 
cancellation of the small divisor 

exp(2i arctan(/iW2/2) + 2i arctan(/iW5/2) — 2i arctan(/iaj_7/2)) — 1 
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which naturally appears in the search for a normal form for the numerical inte- 
grator (2.2) (compare [6, 7]). 

We next consider the midpoint rule defined as 



u 



n+l 



ihHn 



,u 



n+l 



■ ) -thg[ 



n+l 



) 



where g{u) = \u\'^u. Equivalently, this scheme can be written 



,n+l 



R{-ihHo) o [i 



2ih 



ihHn 



))• 



(2.6) 



(2.7) 



We use the same parameters as before. In Figure 3 we observe the preservation of 
the actions for /i = 0.13 and in Figure 4 we see unexpected energy exchanges for 
the step size (2.5). Note that the numerical resonance effect is less pathological 
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Figure 3: Plot of the actions for non-resonant step size (midpoint rule). 



than for the mid-split integrator. This might be due to the regularization effect 
acting on the nonlinearity term (see (2.7)). However the preservation of the 
actions is lost anyway, reflecting the fact that the unperturbed linear integrator 
is the same as for the mid-split integrator, as well as the small divisors to be 
controlled. 



2.2 Resonance-breaker midpoint 

We consider now the same problem, but we set 

CJ2 = 10, = 30 and co^j = 40, 
so that the (physical) resonance relation holds: 

i02 +UJ5 — (jJ-7 = 0. 
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Figure 4: Plot of the actions for resonant step size (midpoint rule). 



The exact solution is plotted in figure 5. We observe energy exchanges between 
the actions. This solution is numerically computed using a splitting method 
under CFL condition with a CFL number close to 1. 

-1.5| 1 1 1 1 1 1 1 1 1 
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Figure 5: Exact solution with resonances. 



In Figure 6, we plot the evolution of the actions of the numerical solution 
given by the mid-split integrator (2.2) with a step size h = 0.0812... such that 

2 arctan(/ia;2/2) + 2 arctan(/iti;5/2) — 2 arctan(/ia-'_7/2) = — . 

We observe the preservation of the actions over long time which reflects the 
control of the small divisor 

I exp(2i arctan(/iW2/2) + 2f arctan(/itiJ5/2) — 2i arctan(/iLt;_7/2)) — 1| ~ 0.485... 
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Figure 6: Undesired actions preservation by the mid-split integrator. 



The same picture is plotted in Figure 7 for the midpoint rule (2.6). We thus see 
that these midpoint-like schemes are both unable to reproduce the qualitative 
behavior of the exact solution in Figure 5. 
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Figure 7: Undesired actions preservation by the midpoint rule. 



Note that the sets of stepsizes h for which the solution looks like in Figure 6 
or 7 for mid-split of midpoint rule integrators are very large. In fact, in this case 
the relation 

j) = ±2arctan(/ia;j^/2) ± . . . ± 2 arctan(/iu;j,./2) = (2.8) 

for a multi-index j = (ji , . . . , jV ) is exceptional, and not implied by a physical 
resonance 

J7(j) = ±wj-, ±...±a;jV = 
unless h is very small so that the linear approximation of arctan(x) is valid. 
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On the opposite, splitting methods will always catch these physical reso- 
nances. 

2.3 Conclusion and summary of the results 

The phenomena described above are somehow generic for any symplectic Runge 
Kutta method applied to Hamiltonian PDEs. Indeed for these methods, the 
unperturbed linear integrator can always be written of the form (2.3) (see for 
instance [10]) and the control of the long-time behavior of the numerical solutions 
thus relies on the control of small divisors 

exp(i^'(/i,j)) - 1 

where ^{h,j) is of the form (2.8). This yields non linear relations between h and 
the frequencies whose control is not implied by the physical small divisors 

To remedy this, a possibility is to use splitting methods (for which ^{h,j) = 
hQ{j)) with suitable truncation in high-frequencies. The goal of this paper is 
to give new results extending the existing ones in [6, 7, 3, 8] by considering 
Hamiltonian PDEs where the operator Hq is replaced by a truncated operator 
such that hujj < K form some K. This yields a truncated Hamiltonian H}^ = 
Hq + P where the nonlinear part is the same as for H = Hq + P. In this paper 
we show: 

• The preservation of the actions for the exact solution of the Hamiltonian 
PDE associated with Hh under a generic non resonance condition on the 
frequencies of Hq. This preservation holds for times of the form e~^h~^ + 

where e is the size of the solution in Sobolev norm. 

• The persistence of this result for the numerical solution obtained by splitting 
methods induced by the decomposition of H^ assuming an explicit bound 
for the cut-off level K. This preservation holds for all h < Hq (without 
additional resonances) and for a number of iterations of the form (1.4). 
The bound on K turns out to be independent of N and r and sharp (see 
Thm. 4.2 for a precise statement). 

We state the results in the Section 4, after recalling in Section 3 the mathe- 
matical framework already developed in [7]. For ease of presentation, we do not 
give the entire details of the proofs, as it uses very technical tools that might 
already be found in [9, 7] and [2]. We prefer to give the key ingredients, and 
stress the specific changes in comparison with [7]. 
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3 Abstract setting 

3.1 Abstract Hamiltonian formalism 

This section essentially describes the formalism used in [6, 7]. We denote by 
J\f = or N'^ (depending on the concrete application) for some d > 1. For 
a = (ai, . . . , ad) G M , we set 



max (l, Oi + • • • + a^) . 



We consider the set of variables (^aj %) G C^'^ x C-'^ equipped with the symplectic 
structure 

iJ^d^aAdr^a. (3.1) 

We define the set ^ = x {±1}. For j = (a, 5) G we define \ j\ = \a\ and we 
denote by j the index (a, —(5). 

We will identify a couple (^, r?) G C-'^ x C-^ with {zj)j^z G via the formula 



ia if 5=1, 

r/a if (5 = -1. 



By a slight abuse of notation, we often write z = (^, ry) to denote such an element. 

For a given real number s > 0, we consider the Hilbert space Vs = ^s(2^, C) 
made of elements z G C"^ such that 



I i|2 
I II ^ 



and equipped with the symplectic form (3.1). 

Let be a an open set of Vg- For a function F of C^(Z//,C), we define its 
gradient by 

'dF\ 



where by definition, we set for j = (a, (5) G x {±1}, 

r dF 



OF 

dzj 



— if 5 = -I. 



Let -^(2;) be a function defined on li. If i7 is smooth enough, we can associate 
with this function the Hamiltonian vector field Xh{z) defined by 

Xh{z) = JVH{z) 



11 



where J is the symplectic operator on Vg induced by the symplectic form (3.1). 
For two functions F and G, the Poisson Bracket is defined as 
, ^ ^^r dF dG dF dG 

We say that z G "Pg is real when z-j = 'Zj for any j G Z. In this case, z = (^,0 
for some i^C^. Further we say that a Hamiltonian function H is real if H[z) 
is real for all real z. 

Definition 3.1 Let s > 0, and letU be a neighborhood of the origin in Vg- We 
denote by Ti^ipi) the space of real Hamiltonian H satisfying 

HeC°°{U,C), and Xh & C°°{U,Vs). 

With a given function H G 7i^{U), we associate the Hamiltonian system 

i = JVH(z) 

which can be written 

dH ^^-^^ 

In this situation, we define the flow (p\j{z) associated with the previous system 
(for times t > depending on z & U). Note that if z = (^, ^) and using the fact 
that H is real, the flow {(,^,rf) = '^^^{z) satisfies for all time where it is defined 
the relation = ff ■, where ^* is solution of the equation 

BTJ - 

ia = -ij^{tO, aGM. (3.3) 
In this situation, introducing the real variables pa and qa such that 

= -^{Pa + iqa) and = -^{Pa - iqa) , 
the system (3.3) is equivalent to the system 

Pa = --k—{q,P) aeM 
oqa 

qa = ■7^{q,p), aeN. 
dpa 

where H{q,p) = H{^,^). 

Note that the flow r* = (/^^ of a real hamiltonian x defines a symplectic map, 
i.e. satisfies for all time t and all point z where it is defined 

{Dy)U{Dy), = J (3.4) 

where denotes the derivative with respect to the initial conditions. 
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3.2 Function spaces 



We describe now the hypothesis needed on the Hamiltonian H. 

Let ^ > 3 be a given integer. We consider j = {ji, ■ ■ ■ ,je) G Z^, and we set 
for all z = ji = {ai,6i) where Oj G and 6i G {il}- We define the 

moment of the multi- index j by 

M{j) = ai6i + ■ ■ ■ + ae6e. (3.5) 

We then define the set of indices with zero moment 

Ie = {j = {ji,---,je) with M{j) = 0}. (3.6) 

For j = (ji, . . . ,jr) G Ir, we define as the third largest integer between 
|ji|, . . . , IjVI- Thenweset S'(j) = - |jv_i | +m(j) where and | J denote 
the largest and the second largest integer between . . . , \jr\. 

In the following, for j = (j'l, . . . , j^) G Ii, we use the notation 



Moreover, for j = (ji,... G le with ji = (a^, 6i) e M x {±1} for z = 1, . . . , 
we set 

3 = Gi,---,Je) with J. = (ai, -(5i), i = !,...,£. 
We recall the following definition from [9]. 
Definition 3.2 Let k > 3, M > and v G [0,+oo), and let 

k 

We say that Q G T^ ^'^ if there exist a constant C depending on M such that 



y£ = 3,...,k, yjele, \Qj\ < C ^l'/ . (3.7) 



SUV 



Note that Q is a real hamiltonian if and only if 

y£ = 3,...,k, \fjeli, Qj = Qj. (3.8) 

We have that T^'" G W for s > + 1/2 (see [9]). The best constant in the 
inequality (3.7) defines a norm \Q\ m.i^ for which '71*^''^ is a Banach space. We 
set 
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Definition 3.3 A function P is in the class T if 

• P is a real hamiltonian and exhibits a zero of order at least 3 at the origin. 

• There exists sq > Q such that for any s > sq, P £ Ti'^iU) for some neigh- 
borhood U of the origin in Vg ■ 

• For all k > 1, there exists > such that the Taylor expansion of degree k 
of P around the origin belongs to T^''^ ■ 

With the previous notations, we consider in the following Hamiltonian functions 
of the form 

H{z) = Ho{z) + P{z) = J2 ^alaiz) + P{z), (3.9) 

aeA/" 

with P £ T and where for all a G 

Ia{z) = ^aVa 

are the actions associated with a G M. We assume that the frequencies G M 
satisfy 

VaGAA, \uJa\<C\ar (3.10) 

for some constants C > and m > 0. The Hamiltonian system (3.2) can hence 
be written 

dP 

= -ii^a^a-i-K—{^,V) a & Af 

(3.11) 

fja = iiOaVa + i-^i^.,'!]) aGj\f. 

3.3 Non resonance condition 

Let j = (ji, . . . ,jr) G Tr, and denote by ji = (cj, 6i) £ M x {±1} for i = 1, . . . , r. 
We set 

= 6iUa, + ■ ■ ■ + SriOa,. (3.12) 

We say that j = {ji, ■ ■ ■ , jr) G depends only of the actions and we write j G Ar 
if j = j. In this situation r is even and we can write 

Vz = l,...r/2, ji = (ai,l), and ji+r/2 = (oi, -1) 

for some Oj G M. Note that in this situation, 

~ ^jl ' ' ' ^jr ~ CaiVai ■ ■ ■ ^a^/2'^o.r/2 

where for all a G A^, 

Ia{z) = CaVa 
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denotes the action associated with the index a. Note that if z satisfies the 
condition zj = 'Zj for all j E Z, then we have Ia{z) = |CaP- For odd r, Ar is the 
empty set. 

We will assume now that the frequencies of the linear operator Hq satisfy the 
following property: 

Hypothesis 3.4 For all r G N, there exist constants 7 and a such that Vj = 

(jl, ■■■,jr) ^ Ar, 

imi > (3.13) 



where we recall that denotes the third largest integer amongst . . . , 

As explained in [6, 7, 2, 9] a large set of Hamiltonian PDEs can be written 
under the previous form and satisfy the non resonance condition (3.13). This 
includes the case of Schrodinger equations on a torus on dimension d and wave 
equations with periodic boundary condition in dimension 1. 

3.4 Numerical integrator 

In the following K is a fixed number and characterizes the high frequency cut-off 
level. As explained in the previous section the integrators we consider are such 
that the frequencies of the linear operator multiplied by h are bounded by K. 
Let xk '■ ]R+ — > ]R+ be the cut-off function 



Xk{x) 



X if X < K 
if x> K 

We define the operator ipKihHo) := exp{ixKihHo)) by the formula 

^-SihuJa if huJa < K 



Vj = (a, (5) e Z, ipK{hHo)j = exp{-i6xK{huJa)) 



1 if huj> K 



We consider split-step integrators of the form 

n"+i=Vi^(/iFo)o^^(^z"), (3.14) 

where ipp is the exact fiow of P. 

The numerical scheme (3.14) corresponds to an exact splitting method with 
time step h applied to the truncated equation associated with the (infinite di- 
mensional) Hamiltonian 

Hh{z)=H^{z)+P(z):= Yl iOa^ar]a + P{z), (3.15) 

a£j\r,\LJa\<Kh-^ 
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where we denote by Hq the truncated hnear operator. The corresponding Hamil- 
tonian system can be written 

1 dP 

L = -i-rXK{huLia)^a-i^{^,v) aeAf 

1 a?" (^-^^^ 

Va = 'i'TXKihuJa)Va + i^{C,V) a^M. 

Note that the effect of xk is a cut-off in uja such that LOa < Kh~^ . 
The following result is easily shown and given here without proof: 

Theorem 3.5 Assume given two solutions z[t) and Zh{t) of (3.11) and (3.16) 

respectively, such that Zh{0) = z{0). Let s,a > and assume that for < t <T 
we have z(t) and Zh{t) in V^^^^^i. Then we have 

vtG(o,r), \\z{t)-zh{t)\\,<ch'' 

where C depends on the Sobolev norms of z{t) in V^j^^j^i for t G (0, T). 

4 Statement of the result and applications 

We first give normal form results both for the truncated equation (3.16) and the 
numerical integrator (3.14). We then give the dynamical consequences of these 
results for the long time behavior of the corresponding solutions. 



4.1 Normal form results 

Theorem 4.1 Assume that P £ T, and that the non resonance condition (3.13) 
is satisfied. Let K and r he given numbers. Then there exist sq and ho such that 
for all s > sq there exist Us and Vs two neighborhoods of the origin in Vs such 
that for all h < ho there exists Tg :Vs ^ Us « canonical transformation which is 
the restriction to Us of t = t^q and which put the Hamiltonian of eqn. (3.15) 
under normal form 

HhOT = H^ + Z + R 

where Hq is the Hamiltonian defined in (3.15) and where 

(i) Z is a real hamiltonian, polynomial of order r in z with terms that either 
depends only on the actions or contain at least three components Zj with 
\j\ > Kh~^. As a consequence we have 

yzGVs, J]|a|2^|{I,,Z}(z)|<C7/i^(||z||^ + ||.|Q (4.1) 

aGA/" 

where C depends on r, s and K. 
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(ii) R G WiVs^'C) is a real hamiltonian such that for z € V^, we have 

\\Xr{z)1 < C\\z\\l (4.2) 

where C depends on r, s and K . 
(in) T is close to the identity in the sense where for all z £ Vg we have 

||r(z)-z|| <C||z||; (4.3) 



and for all z £ 



\T-\z)-zl<C\\z\\l (4.4) 



where C depends on r, s and K . 

As explained in the proof, this result is a mixed between the truncature system- 
atically made in [2] and the global result stated in [9]. The dependancy on h in 
the estimates reflects the control of the non resonance conditions associated with 
the truncated linear operator appearing in H^. 

Theorem 4.2 Assume that P £ T, and that the non resonance condition (3.13) 
is satisfied. Let r be a given number. Let K be such that 

K<\. (4.5) 

then there exist sq and Hq such that for all s > sq there exist Us and Vs two 
neighborhoods of the origin in Vs such that for all h < ho there exists a canonical 
transformation Ts : Vs ^ Us which is the restriction to Ug of t = Ts^^ such that 

r"^ o iljKihHo) oip),oT = ipKihRo) o 9 (4.6) 

where 6 is the solution at time X = 1 of a non- autonomous hamiltonian hZ(\) + 
R{X) with 

(i) Z{\) a real hamiltonian depending smoothly on X, polynomial of order r in 
z with terms that either depend only on the actions or contain at least three 
components zj with \ j\ > (r — 2)~^ Kh~^ . 

(ii) R{X) £ 7i^{Vs,'C) a real hamiltonian depending smoothly on X £ (0, 1), and 
satisfying (4.2) uniformly in X. 

(Hi) T close to the identity in the sense where it satisfies (4.3) and (4.4). 

As a consequence, there exist a constant C depending on r and s such that 

^z£Vs, ^|a|2^|4(e(z))-/,(z)| <C(/i^||z||' + ||z||;+'). (4.7) 

Remark 4.3 As will appear clearly in the proof, the bound (4.5) can be refined 
to 11/2 — 6 with 6 > for general situations. If the bound tt/2 is not satisfied, 
we can construct a system such that numerical resonances between h and the 
frequency vector to appear. 
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4.2 Dynamical consequences 

We now give the main outcome of the previous theorems: The first concerns the 
exact solution of (3.16) and the second the long time behavior of the numerical 
solution associated with splitting methods applied to this equation. 

Theorem 4.4 Assume that P G T and Hq satisfies the condition (3.13). Let 
r,N £ N* be fixed. Then there exist constants sq and Hq depending on r and N 
such that for all s > so, there exist a constant eq depending on r, N , K and s 
such that the following holds: For all e < Eq, h < ho and for all £ Vg real 
such that \\z^\\^ < e, then the solution Zh{t) of (3.16) with z/i(0) = z^ satisfies 

\\zh{t)\\^<2e for i < + (4.8) 

and 



Y,H''Va{zH{t))-Ia{zHm\<e"^ for t < ^ + (4.9) 

for some constant c depending on s, r and N . 

Note that Eqn. (3.16) is a infinite dimensional PDE. The only frequency 
cut-off is made in the linear operator. The difference with the classical results 
[2, 9] is the dependence of the cut-off parameter in the bound in time. For 
fully discretized systems obtained by pseudo spectral methods, the same result 
holds with constant independent of the spatial discretization parameter (a priori 
independent of h and K). We do not give the proof here and refer to [6] for the 
description of fully discretized systems. 

Theorem 4.5 Assume that P £ T and Hq satisfies the condition (3.13). Let 
r,N £ N* be fixed, then there exist constants sq and ho depending on r and N 
such that for all s > sq, there exists a constant eq depending on r, N and s such 
that the following holds: For all e < £q, h < ho and for all z^ £ Vg real and 
\\z^\\^ < e if we define 

z^ = {i,K{hHQ)oip%y{z^), (4.10) 
where the frequency cut-off is such that 

K<-. (4.11) 
then we have still real, and moreover 

Wz'^W <2e for n < ^ + (4.12) 
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and 

Y,\a\'Va{zn-Iaiz')\<e'/' for n < + (4.13) 
for some constant c depending on r, N and s. 

Proof of Theorem 4.4. Let = t~^{z^) which is well defined provided £q is 
small enough so that ^Ug- Using (4.4) we have 

SO that we can assume that \\y^\\^ < e/2. We then define yh{t) the solution of 
the Hamiltonian system associated with the Hamiltonian Hq + Z + R given in 
Theorem (4.1). We set M{t) = \\yh{t)\\l = J^sf'^jiVhit))- We have 

s 

Assume that sq > N and using (4.2) and (4.1), we see that as long as y/i(t) £ 
we have 



dt 



< 



C(/i~AA(t)t+AA(t)^) 



< C(/i^e3 + e'-+i^ 



as we can always assume that ho < 1 and that is contained is the ball of radius 
1 inVs- We know that M{0) < and we can assume that eq is sufficiently small 
in such way that the ball in Vg centered at the origin and of radius 3e/2 is 
included in V^. Now as long as ■\/M{t) < 3e/2 we have 

for some constant C. This implies 

\Af{t) - Af{0)\ < Ct{h^e^ + e^'+i). (4.14) 

Hence there exists a constant c such that as long as 

t < c{h-^ e-^ + e'-^) 

we have ^^M{t) < 3e/2 and hence yh{t) G V^. But this implies that Zh{t) = 
T{yh{t)). Using (4.3) we then easily see that (4.8) is satisfied. The proof of (4.9) 
is similar (see [2, 9]). ■ 

Proof of Theorem 4.5. The proof follows the same lines as above but on a 
discrete level. Setting y^ = t~^{z^) we have \\y^\\ ^ < s/2. Now if we define by 
induction 

y''^^ = i,K{hHo)oe{y^) 
and if we define Mn = \\y^\\ ^ then estimate (4.7) shows that 

\Mn+l-Mn\ < C(/i^AAJ +J\f^), 
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provided sq > A^, which is the discrete version of (4.14). We then easily conclude 
upon using the same arguments. ■ 



5 Proof of the normal form results 
5.1 The continuous case 

We start now the proof of Theorem 4.1. 

The strategy follows lines of [2, 9]: we search by induction a transformation 
eliminating the polynomial term of order i until i = r. The transformation r is 
then defined by the composition of all these transformations. At each step, the 
transformation is constructed as the flow at time 1 of a Hamiltonian system with 
unknown Hamiltonian x- Hence we are lead to solve recursively the homological 
equations 

{H^,X} + Z = G (5.1) 

where 

G = ^ GjZj G Tg"^'"^ 

is a polynomial of order i depending on the term constructed in the previous steps 
and satisfying estimates of the form (3.7) for some z^i > 0, and Z an unknown 
term in normal form. 
Setting 

OJa if ^ Kh" 

(5.2) 

if > Kh-^ 



we see that the equation (5.1) can be written in terms of the coefficients Xj^ Zj 
and Gj as 

^^{3)Xi + Zj = Gj (5.3) 

where ^^{j) is defined as (3.12) with respect to uj^. Following [9] in the case 
where no frequency cut-off is made (i.e. ^^{j) = ^{j)), the condition (3.13) 
ensures that the system (5.3) can be solved by putting the terms depending on 
the actions in Zj and solving the rest to construct Xj by inverting r2(j). It is 
then clear that Xj belongs to some T,°°'"^ for some 1^2 which ensures the control 
of the regularity of the transformation. 

In our situation, it is clear that does not fulfill the condition (3.13): if 

for instance (ji,j2)j3) is a multi-index with all components greater than Kh~^ 
in modulus, then ^^{j) is equal to zero. 

On the other hand, let j = (ji, . . . £ Ie\Ai a multi-index with at most 
two indices greater than Kh~^. We can always assume that these two big indices 
are ji and j2 with ji > j2 and that = js- 
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If both are greater than Kh ^ we have in fact 

\^\j)\ = \^{h....,u)\> 



upon usmg (3.13) unless (js, ■ ■ ■ jji) G Ai-2- But in this last situation, the 
condition Al(jf) = implies that ji = j2, i.e. j £ At. 

If only one is greater than Kh~^ then we have with similar notations 

\^\j)\ = m2....,n)\> 



thanks to (3.13) unless (js, ■ ■ ■ ,je) G Ai^i. If i is even this is impossible. 
If i is odd, then the zero moment condition = implies that ji = 

which is a contradiction. 

This shows that (3.13) holds for except for j G Ae or for j such that 

at least three indices are greater that Kh~^ in modulus. Hence we solve (5.3) 
by defining Zj = Gj and Xj = when j G Ai or when j contains at least three 
indices are greater that Kh~^ in modulus; while Zj = and Xj = ^^{j)~^Gj in 
the other cases, i.e. when j ^ Ai and j contains at most two indices are greater 
that Kh~^ in modulus. 

Now the condition (3.13) ensures that Xj ^ T^''^'^ for some 1^2 • The conclusion 
now follows [9] except the derivation of equation (4.1) which is a consequence of 
Lemma 4.11 in [2]. 

5.2 Splitting methods 

We prove now Theorem 4.2. 

We follow now the methodology developed in [6, 7]. We embed the splitting 
method ^'x(/ii7o) o into the family 

(0,1) 3 \ ^4,K{hH^)o^lp 

and we seek as a transformation r = V^(^x^ associated with a non-autonomous 
real hamiltonian xW depending smoothly on A and such that for all A G (0, 1) 

V A G [0, 1] i^KihHo) o iflp o (^^(^) = (^^(^) o Vi^(/iFo) o (5.4) 

where Z{\) is a Hamiltonian in normal form in the sense of Theorem 4.2 and 
R{\) a real Hamiltonian possessing a zero of order r + l. Deriving this expression 
in A, we find the equation (compare eqn. (5.18) in [7]): 

VAg[0,1] xWoi^KihHo)-xWoVhP = hP-{hZ{\) + Ri\))o^-^^^. (5.5) 
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As in [7], we see that the solution of this equation rehes on the solvabihty of a 
discrete Homological equation of the form 

V j G T,, (e^'^^'^^') - l)xj + hZj = hGj 

where now ^^{j) is defined as in (3.12) with respect to lo^ defined in (5.2). 
Assume that Gj G j'^''^^ for some I'l. 

Let j £ ^ Ag. As before, we assume that \ji\ > \j2\ > IjsI > • • • • We recall 
that for |x| < vr, we have 

|e'^ - 1| > -\x\. 
vr 

Assume that /i(j) = js < (l - 2)~^Kh~^. Then we have 

\hn{h,...,je)\<K<l- 

Now we have three possibilities: 

• 1^2! > Kh~^ . In this situation we have uj^^ = a;^^ = 0- 

\Q.\3)\ = \hn{h,...,u)\<^ 

and hence 

2.^^,. hi* 



for some constant 7* using (3.13) and unless {j3,---,je) £ A.e-2- But in 
this last case, thanks to the the zero moment condition, ji = j2 and thus 

> Kh~^ and |j2| < Kh~^. Now we have 



and hence 



2tt 

\n\j)\ = \hn{j2,...,h)\<2K< — 



.^'^'(^^-i\>hn\j2,...,n)>^ 



unless {j3,---,ji) G Ai-2 which is impossible because the zero moment 
condition would be violated. 

\ji\ < Kh~^. In this case we have 

\^\3)\ = < 3J^ < vr 

and hence we have 



jhnHj) _ i| > •) > Jfr 



using (3.13). 



22 



So far we have proven the foUowing: For all j such that fi{j) < {i — 2) ^Kh ^ 
and j ^ Ae, we have 

UhQ^j) _ ^1 > hj* 

~ Kj)" 

for some constant 7* depending on the constant 7 in (3.13) (for different i). 

Now we see that we can solve the homological equation for those indices, and 
put the remainder terms in Z, which entails into the definition of normal form 
terms in Theorem (4.4). It is then clear that xW ^ T^''^^ for some 1^2, and the 
rest of the proof can be adapted from [9] . 

Finally, estimate (4.7) is a direct consequence of Lemma 4.11 in [2]. This 
concludes the proof. 
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